Thermodynamics of nano-cluster phases: a unifying theory 
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We propose a unifying, analytical theory accounting for the self-organization of colloidal systems in 
nano- or micro-cluster phases. We predict the distribution of cluter sizes with respect to interaction 
parameters and colloid concentration. In particular, we anticipate a proportionality regime where 
the mean cluster size grows proportionally to the concentration, as observed in several experiments. 
We emphasize the interest of a predictive theory in soft matter, nano-technologies and biophysics. 



There is an increasing interest in colloid science and 
emerging nano-technologies for systems exhibiting self- 
assembled, nano-structured spatial patterns at equilib- 
rium (pseudo-periodic structures in 2D or 3D, such as 
stripes, ripples, bubbles, lamellae, tubes or clusters in a 
large variety of systems) P, 01 • A competition between a 
short-range attractive interaction {e.g. a depletion force) 
which favors condensation, and a weaker, longer-range 
repulsion {e.g. electrostatic) which prevents a complete 
phase separation, leads to equilibrium structures with 
a characteristic length-scale (the pattern typical size). 
Here we focus on cluster phases, a particular type of space 
patterns, consisting of aggregates of sizes ranging from a 
few to hundreds of particles. They are the colloid analog 
of micelles Q and occur for a variety of physical sys- 
tems {e.g. colloids, star polymers, proteins) and interac- 
tions [11, i, i, S, 0, S, i, [13, El, [H , whenever the strength 
of the interactions is of the order of the thermal energy 
/jfiT. Such patterns have also been experimen tally ob- 
served or simulated in two dimensions [3, [13, 12, [13] ■ 
Recently, it has been proposed by one of us that proteins 
embedded in cell membranes can also be found in two- 
dimensional cluster phases [l4|. Indeed, they experience 
repulsive and attractive forces mediated by the lip idic 
membrane, the range of which is nanometric [l5, 13, l3|- 
Nowadays, such cluster structures can be observed in cell 
plasma membranes with the help of advanced microscopy 
techniques 13, i3|- 

It is of great interest in this context to be able to pre- 
dict a priori for which regimes of parameters one can ex- 
pect a cluster phase to exist, and to anticipate the mean 
cluster size. However, apart from generic phenomeno- 
logical arguments [2I, \a{ , the thermodynamics of cluster 
phases has only been tackled in specific frameworks and 
dimensions, ranging from Van der Walls approximative 
schemes (3] to electrostatic models taking explicitly into 
account counter-ions [3, [131 , or estimations of structure 
factors U, 2l[. By contrast, we propose here a general 
approach based on simple thermodynamical principles, 
requiring no specific description of microscopic interac- 
tions, but simply the existence of some generic features 
of these interactions, and applicable in two and three di- 



mensions. The analytical methods appeal to ideas from 
elementary micellisation or nucleation theories [3, [23 |. 
They are also reminiscent of Ref . [23| , even though this 
study was limited to steric repulsion between membrane 
proteins (see also [l3|)- Up to basic prerequisites, we 
demonstrate that the cluster phase exists above a critical 
particle volume fraction, 0"^ (Fig. [1]). Then a gas phase 
coexists with large clusters. In addition, in several cir- 
cumstances, the mean aggregation number (fc) has been 
measured experimentally in function of the volume frac- 
tion (j); (k) is extracted either from direct enumeration 
in confocal microscopy [ll 13 or from structure factors in 
diffraction experiments -"-^ found that for a wide 

regime of concentrations, (fc) grows proportionally to (f): 



{k) ~ (/./^^ 



(1) 



Our analytical treatment provides a straightforward ex- 
planation for this proportionality regime (Fig. [1]). 
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FIG. 1: Exact mean aggregation number (fc) iis colloid volume 
fraction for different sets of parameters, in log-log coordi- 
nates, calculated within our framework (Eqs. (|6l7|l 'l. Here 
d = 2, a = 3/2, 7 = 40 and / and a are indicated in the 
legend (see definitions in the text and Eq. (I14p ). The dotted 
line has slope 1, for comparison. For (f) < (f)'^ , the systems 
is in the gas phase and (fc) ~ 1. The proportionality regime 
where (fc) cx appears clearly for (j)'^ < <j) ^ '/'sup, followed by 
a saturation for tj) > "i^sup where large multimers dominate. 
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The cluster size distribution is investigated in a statis- 
tical mechanics formalism. A useful introduction to this 
issue can be found in the Appendix of Ref. [2§| . We start 
from the canonical partition function for a system of N 
interacting particles in a volume V of dimensionality d, 



Z = 



A 



-dN 



dri . . . drjv e 



-Un 



(2) 



where Un is total interaction energy of the iV particles; 
here and in the sequel, all energies are in units of ksT. 
Usually, A is chosen equal to the de Broglie thermal wave- 
length A = 2Trh^ /uikBT. However, for sake of conve- 
nience, we shall choose it to be the particle diameter. 
Suppose now that one is able to determine unambigu- 
ously regions Vk of V that partition the A'' particles into 
clusters: A^i monomers, N2 dimers, etc. . . , each k-mer in 
a distinct region Vk , so that N — J2 kNk ■ Then the in- 
tegral ([2]) can be written as a sum of integrals on the Vk , 
because interactions between the regions are negligible. 
Simple algebra 22] leads to 



d^-F(k) 



(3) 



{Nk} k 

where Fll) — and 



Fik) = - In 



fc! 



dri . . .drfe_i e 



(4) 



for fc > 1. It is the free energy of a fc-cluster. 

The particle organization is described by the mean vol- 
ume fraction of fc-clusters which derives from Z, 



Ck 



{Nk)A' 



V 



cie 



-G{k,fi) 



(5) 



where we have introduced the chemical potential /i = 
Inci and the grand potential of a fc-cluster C?(fc,/i) = 
F{k) — (fc — The value of /i, or equivalently of the 
monomer fraction ci , is fixed by the constraint [s, [2^ 



^fcCfe 

k=l 



(6) 



We also define the total volume fraction of clusters (in- 
cluding monomers) M, the total fraction of multimers M 
and the mean cluster aggregation number (fc) , 



00 
fc=i 







00 

Cfc , A/ = ^ Ck and (fc) 

fe=2 



(7) 



The interaction energy of the fc particles of a cluster, 
Uk, is the sum of two contributions: Uq due to the short 
range interactions (hard core repulsion at contact and 
short range attraction) and U-c due to the longer range 
repulsions. The free energy can be split in two parts. 



F{k) = fo(fc) + F^{k); Fq is the free energy given by 
(4) in the absence of repulsion {U = Uq) and Fr(fc) = 
-ln(e-^')o, with (...)„ = e^"W / dri...drfc_i (...)e-^o, 
is the specific contribution arising from the repulsion. 

If particle repulsion is switched-off, the cluster free en- 
ergy reduces to that of a spherical droplet of simple liq- 
uid. It is usually written within a good approximation 
as the sum of a bulk and surface energy, for fc 3> 1: 



^^o(fc) = -/o(fc-l) + 7(fc-l)" 



(8) 



The free energy per particle — /o < accounts for 
the mean energy eo ~ —eav/'2, where v is the typ- 
ical number of neighbors of a particle in the cluster 
and Ea the strength of the short range attraction poten- 
tial between two particles, and an entropic contribution 
— fcs [ln(w/A~^) -I- 1] {vf is the free volume per particle 
accessible to the particles inside the clusters). For inter- 
mediate, finite, values of fc, a positive surface correction 
must be added to the previous bulk contribution. It takes 
into account the fact that the surface particles have typ- 
ically twice fewer neighbors than the bulk ones and a 
larger free volume. 

Now we consider that a weak repulsive pair potential 
is also acting on the particles. 



C^r(fc) = i5]^;(ri-rj). 



(9) 



If we assume that those weak interactions do not modify 
the cluster structure and that the particle density, p = 
kA'^/Vk « 1, is homogeneous in the cluster, we can make 
the two successive mean-field approximations: 

F-ik) (t/,.)o =^ (fc - 1)^ / dridr2 v{r^ - r^). 

^Vk JVkXVk 

. 

As an example, we consider the following potential of 
intermediate range A > A, 



V{r) = Er ( - 

1 r 



(2-a) 



-r/A 



(11) 



For r < A it behaves like a long range potential if a > 1. 
But at larger distances, the potential rapidly vanishes. 
For particles undergoing screened Coulomb repulsion, 
A is the Debye screening length, a = 5/3 in 3D and, 
a = 3/2 if the particles are confined in 2D, like charged 
proteins in a membrane. Using the potential (jlip . by 
estimating the integral in the repulsive part of the free 
energy (fTO| . one gets Fr{l) = and. 



F,{k)^f,ik- 

Fr{k) - 



-l) + a{k- 

/r(fc-l), 



1)", fc<C pA^A"^, (12) 
fc » pA^A^ (13) 



with cr - £i.(pl/'iA/A)(2-")'i. 
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First we consider that the repulsion is of infinite range 
(A ^ oo) so that p2| holds for any k. According to Eqs. 
(|8ll2p the grand potential reads 

G(fc,M) = -{f + fi)(k-l)+j{k-l)^+<j{k-ir, (14) 

with the bulk free energy / = /o — /r- The existence of 
clusters even at significantly low concentrations {(j)'^ <^ 1) 
requires / 3> 1 to overcome translational entropy (see 
proof below and Fig[T]). We also assume that the repul- 
sion is weak: e,- <C 1 so that a < 1 and f ^ a. 
Finally, we shall see below that the existence of the pro- 
portionality regime requires 7 3> cr. The global shape of 
G{k) - and thus of Ck, Eq. ([5|) ~ is very sensitive to the 
value of /i, which is itself fixed by the total concentration 
of particles (Eq. [H])); Eqs. (|5l6p show that increases 
monotonously with (p. The particle concentration thus 
controls directly the distribution shape. 

At low particle concentration (f), = Inci < ln</) 
takes a large negative value and then, G{k) increases 
monotonously with k. It follows from ([5]) that the clus- 
ter size distribution Ck is maximal at k = 1 and decreases 
exponentially with k. Since the behavior of G{k) is dom- 
inated by ~fik, the width of the distribution Ck goes as 
The volume contains mainly single particles and 
very few transient and small clusters formed by thermal 
fluctuations. The particles form a gas phase. The mean 
cluster size, (fc) ~ 1, is slightly above 1. 

Above a critical concentration {(p > (j>'^), the chemical 
potential /i > //"^ is such that G{k) has a local minimum 
at A:* > 1. In this regime, the distribution Ch is bimodal 



with two maxima at A; = 1 and k — k* 1J|. The par- 
ticles are partitioned between a gas of monomers and 
stable clusters of aggregation number fluctuating about 
k* . Such a conflguration is often called a cluster phase. 
Since the distribution is peaked at k* , we can estimate 
with a good accuracy the multimer concentration 



M ~ cie 



(15) 



and we write M = ci + M, (p ~ ci + k*M . 

At the critical point fi = jj,'^, G{pf , k) has an inflexion 
point at k = k* = k'^. It follows that /x'^ and k'^ satisfy 
dkG{k'' , n") = dlGik^.ii") = and then read. 



r 1 = 



1 d- 1 7 
a(a — 1) a 



(16) 
(17) 



where A > is a long prefactor, function of a and d. The 
condition 7 3> cr ensures that fc"^ ^ 1. The requirement 
(j)'^ implies ^l"^ — \nc\ < ln0^ ^ — 1. Eq. (fTTl) implies 
/ > — /Lt'^, which justifies the condition / ^ 1 above. In 
addition, G(fc°, n") » 1 thus M" < and 



In the cluster phase close to the critical point, k* ~ k'^ 
and, since G{k'^, fi'^) ^ 1, the concentration of monomers 
exceeds by far that of multimers, M <C ci (fTS]) . As we 
go deeper into the cluster phase by increasing </> or /j,, 
the typical cluster size k* shifts to larger values and M 
raises as the energy well G{k* , fi) deepens. In order to get 
quantitative insights on the behaviors of the cluster phase 
above the critical point, among which the demonstration 
of the proportionality law ([T]), we perform a systematic 
expansion in terms of the small parameter (fc* — k'^)/k'^ > 
0. Since those calculations involve very long prefactors in 
the general case, we treat here the particular case d = 2 
and a = 3/2 [3, [IJ, |23|; for other values of d and a, the 
calculation can easily be performed following the same 
route. In this case, we have fc^ = 1 + 7/3(T ~ 7/3cr ^ 1, 

Ai" = -/ + and G(fc^Ai") = i^f ^ » 1- First 

the deviation of the chemical potential from its critical 
value is obtained by doing the expansion of the equation 
satisfied by k* at a given /i, dkG{k*,fi) = 0, in terms of 
^JL~ ^jf and k* — k'^ up to the first significant order. Thus 



^ %/3 , fk* 

- - ^\rw 



(19) 



Next, using we can obtain the increase of the cluster 
free energy, AG = G{k*,n) - G{k'^,^''), as 



AG~-^G(r,/i=) ' 



(20) 



It allows us to express the increase of the monomer and 
multimer concentrations from their values at the critical 
point cl = e^"" and M" = cJe-^('=''^') 

ci-c^ c.c^(a.-A^=), (21) 
M - M" ~ M^i^i- ^f + e-^^ -I). (22) 



According to Eqs. (|19ll22p . we finally obtain the increase 
of the total cluster concentration M = ci -\- M: 

M -M" = Cl - Cl + M - M", (23) 
~ M'=(M-/i^) + M'=(e-^^-l), (24) 

as well as the increase of ~ ci -I- k*M: 

(j) ~ Cl - + k'={M - M'^) + M'={k* - k'=){25) 



-AG 



1 



k* - k' 



(26) 



(18) 



We have expanded G in powers oi k — k'^ at the lowest 
significant order. One can prove that these expansions 
are relevant while k — k'^ < k'^. In contrast, expansions 
of e~'~^ become rapidly erroneous when k ^ k'^ = 0{k'^). 
Thus we keep the exponentials in the expansions. 

The inspection of Eqs. ()20|24|26|) reveals that when 
k* — k'^ grows, the total fraction of colloids, (j), can in- 
crease considerably without the total fraction of clusters. 
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M, varying significantly. It follows that (k) = (j)/M 
grows linearly with (p within a very good approxima- 
tion, which proves the proportionality law More 
precisely, M — M'^ remains of the order of M'^ provided 
that (i) ti- ^J,'' < I, i.e. k* - < [d,/ ^/i^f^k^- 
imd (ii) M^e-^^ < M", i.e. Gik^.n") + AG > or 
k* — k'^ < 2\/2/3fc^. We denote by kgnp the ensuing lim- 

1/2 ^1 

3 



iting value: k. 



sup 



fc'^ = inf 



\V3^ J ' 3 



fc^ For 



reasonable parameter values a < 0.5 and 7 < 100, such 
as in Fig. [2 the numerical prefactor above is close to 1 or 
larger. We simplify below this condition to fcgup = 2fc^. 

By contrast to M, <f> grows rapidly with ^ because of 
the prefactor k'^ in the second term of Eq. (pS)) . When 
k — fcsup, one gets (/)sup — ^^0^ by Eq. (^5]) . Thus 
'/'sup/0'^ = fcsup/2 — (fc), because at this (j), about one half 
of the clusters are monomers and the other half are mul- 
timers of size fcsup- Therefore the proportionality regime 
remains valid up to </> = (/fsup, as illustrated in Fig. [T] 

To sum up, there are 3 regimes: ii (p < (j)'^ di — 
e'^ , the system essentially contains monomers; if > 
(/)sup, far deep in the cluster phase, where the fraction of 
clusters exceeds the concentration of single particle, M > 
ci, the average aggregation number is obviously (fc) ~ 
k* . It varies hardly with (j). Finally, the proportionality 
regime is the range of concentrations (jf < (j) < 4>sup 
where (fc) ~ 4'/4''^ — 4': in other words, owing to 
Eq. (O, where the cluster volume fraction M varies very 
slowly. The proportionality regime remains valid while 
the distribution Ck is truly bimodal, that is to say the 
fraction of multimers remains comparable to the fraction 
of monomers j24| . This regime is visible in Fig. [T] An 
estimate of </>sup is (/>sup — — ^'f^'^ ^ 'Z'^- The 

proportionality indeed covers about two decades in Fig.[T] 

It is also worth discussing the role of small multi- 
mers that we have neglected so far. In our framework, 
monomers are far more numerous than dimers (or small 
multimers): C2/C1 = exp(/x'^ — fj) exp(7 — f + a — ^) — 
exp(^'^ — /i) exp(7 + a — ^30-7) > exp(/x'^ — fi) exp(7/4). 
Thus C2/C1 3> 1 if /i ~ fj.'^, because 7 S> 1. This impor- 
tant point justifies the use of the grand potential p4| for 
small clusters. Indeed, the form of was a priori valid 
for large k only, where the definition of a surface tension 7 
and of a "bulk" energy /o is meaningful. A correct mod- 
eling of dimers would involve the true binding free energy, 
Fh{2), of a single bond between monomers. With Fi,{2) 
of a few ksT, one also gets that C2/C1 = el^''(^)l+'' < 1, 
because /i <C — 1. The same arguments hold for Cfe/ci 
when 1 < k <^ k'^. Small multimers are negligible and 
their exact modeling is irrelevant in our formalism. 

Finally, we now explore the case of physical interest 
A finite. Eq. ([13]) shows that at large k, F recovers a 
simple spherical droplet form, as in Eq. ([8|). If / = /o — 
fr > 0, then the equilibrium configuration at the large N 
limit is a single large cluster coexisting with gas, because 
this configuration saves the surface free-energy cost (if 



/ < 0, the condensation is never favorable). Inspecting 
Eqs. (|12m4p . if there exist k - X'^/A'^, k > k* where G is 
maximum and such that G{k,ii) — G{k* , fi) 3> fcsT, the 
cluster phase becomes metastable. 

Thus we have shown that the different physical quan- 
tities of interset can be predicted by our analytical ap- 
proach. Starting from the microscopic interaction po- 
tential, we predict a range of concentration where (i) a 
gas of monomers coexists with large clusters, and (ii) the 
mean aggregation number is proportional to the colloid 
concentration, as observed in experiments. We believe 
that, playing with the effective interaction parameters, in 
particular the colloid electrostatic repulsion, the Debye 
screening length and the attractive part due to depletion 
forces, our predictions can be tested experimentally by 
confocal microscopy where cluster statistics can be ob- 
tained with a very good accuracy [l|. 
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